Simulation of Early-stage Clustering in Ternary Metal Alloys Using the Phase Field Crystal Method 
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Phase field crystal methodology is applied, for the first time, to study the effect of alloy composition on 
the clustering behavior of a quenched/aged supersaturated ternary Al alloy system. An analysis of the work 
of formation is built upon the methodology developed in Fallah et al, to describe the dislocation-mediated 
formation mechanisms of early clusters in binary alloys [Phys. Rev. B., DOP 10. 1 103/PhysRevB. 00.004100]. 
Consistent with the experiments, we demonstrate that the addition of Mg to an Al-l.lCu alloy increases the 
nucleation rate of clusters in the quenched/aged state by increasing the effective driving force for nucleation, 
enhancing the dislocation stress relaxation and decreasing the surface energy associated with the Cu-rich Cu- 
Mg co-clusters. Furthermore, we show that it is thermodynamically favourable for small sub-critical clusters to 
have higher affinity for Mg than larger overcritical Cu-rich clusters, particularly depicting a two-stage clustering 
phenomenon. 



I. INTRODUCTION 

The earliest stage of structural decomposition during 
quench/ageing of supersaturated solid solutions, referred to 
as solute clustering, is a crucial step for establishment of their 
final microstructure. Solute clustering controls the mechani- 
cal properties of alloys through the dispersion pattern of small 
coherent lattice aggregates, namely clusters and/or GP zones. 
This phenomenon, also known as early-stage age hardening, 
is highly influenced by the chemical composition of the alloy, 
especially in multi-component systems. The physical mecha- 
nisms of clustering have been poorly understood by investiga- 
tors exploring the effect of adding different elements to binary 
and ternary systems, as the atomistic behavior of the cluster- 
ing phenomenon are often challenging to model or fully char- 
acterize at the atomic scale. 

Understanding the atomistic mechanisms of the clustering 
phenomenon in multi-component alloys is of crucial impor- 
tance to efficiently design age-hardening processes for desired 
properties, and to accurately interpret the experimental obser- 
vations and measurements. A systematic study of the clus- 
tering mechanisms precludes traditional atomistic methods, 
such as molecular dynamics (MD) and thermodynamic Monte 
Carlo (MC) simulations, which cannot operate on the diffu- 
sional time scales controlling clustering and related solid state 



transformations. Dynamical calculations with classical den- 
sity functional theory (CDFT) also do not apply since they too 
operate on too small a time scale to be relevant to diffusion- 
controlled phase transformations processes [ 1 1. 

Recently, a formalism coined the phase field crystal (PFC) 
methodology 12^61 has emerged that contains many of the 
salient fundamental principles of CDFT but which is suit- 
ably simplified so as to render calculations of microstruc- 
ture kinetics with atomic-scale effects tractable on diffusive 
time scales. The atomic density field in the PFC formalism 
is coarse-grained in time [ 1 1 and does not have sharp peaks 
in solid phases allowing lower spatial resolutions. Numer- 
ous studies have reviewed and demonstrated the physics of 
the PFC methodology and its usefulness in describing a range 
of non-equilibrium microstructure phenomena, from solidifi- 
cation [7, 8 1 and grain boundary kinetics J9] [10] to cluster- 
ing ifTTIl and phase patterning due to atomic misfit strains lfl2l . 
The most recent PFC formalism developed by Greenwood et 
al. (5] [6) employs correlation kernels in the free energy that 
stabilize various crystal symmetries, and coexistence between 
them. More recently the approach was extended to binary al- 
loys, represented by the dynamics of density and a concentra- 
tion field Q3) . 

In our most recent investigation ifTTl with the binary PFC 
model of ref. fl3l . we systematically elucidated a com- 
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plete free energy path for early-stage clustering mediated by 
quenched-in dislocations in the bulk crystal, a mechanism in- 
ferred initially from previous experimental findings in binary 
alloys lfT4l IT31 . In particular, we showed that the energy bar- 
rier for formation of stable clusters can be lowered or even 
completely removed locally in the bulk matrix in the presence 
of an assembly of quenched-in dislocations. Here, we extend 
our energy analysis to ternary systems. The Al-Cu-Mg system 
is chosen since it has been vigorously studied for the evolu- 
tion of clusters [16-21] and shown to exhibit enhanced clus- 
tering and age hardening by Mg alloying |fl6l [181 . Moreover, 
the dominant effect of elasticity and, more particularly, the 
role of quenched-in dislocations has been observed during the 
early-stage decomposition of these alloys in quenched/aged 
state 11201 . 

In this paper we use a newly developed ternary phase 
field crystal model to study solute clustering phenomenon in 
ternary Al-Cu-Mg alloys. The details of a multi-component 
PFC model are presented in a separate paper. In this work, 
we focus on the effect of addition of a ternary species, i.e., 
Mg, on the clustering behavior of these alloys. The remain- 
der of this paper is organized as follows. Section [IT] begins 
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with an introduction of the ternary PFC model. Section 
then demonstrates the model's equilibrium properties, explic 
itly focusing on the Al-Cu-Mg system. Section IV then dis 



cusses new simulations showing the microstructural and com- 
positional evolution of clusters. Section[V]details the analysis 
of system energetics in terms of cluster composition and work 
of formation during evolution. Where appropriate throughout 
Sections IV and [V] results are compared with experimental 
data in the literature. 



II. TERNARY PFC MODEL 

This section reviews the main features of a simplified three- 
component free energy functional from which the dynamics 
of a phase field crystal density field and two impurity concen- 
tration fields are modelled. The starting point of our model 
is a multi-component analogue of the classical density func- 
tional theory (CDFT) of freezing introduced by Ramakrishan 
and Yussouff (22). Details of the derivation of the multi- 
component model are presented in a separate publication 11231 . 
Only the details relevant to the ternary model used in this 
study are reproduced here. 



A. Simplified ternary PFC Free Energy 

The free energy functional of a three-component system 
can be described by two contributions, ideal and excess en- 
ergy, each as a function of three density fields (i.e., pa, ps 
and pc)- From CDFT, we can write the following energy 
functional: 



AF 

k B Tp' 



'^y= J drf = Jdr {AF ld + AF ex }, 



where AF^ are, respectively, the dimensionless the ideal en- 
ergy, AF ex and excess energy, where k B the Boltzmann con- 
stant, T the temperature, p° the average density (defined be- 
low) and V the volume of the unit cell. The ideal free energy 
of the mixture is given by 
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where pi (with i = A,B, C) is the density of component i, 
5p = pi — p° and p° the reference density of component i, 
taken to be that of liquid at solid-liquid coexistence. 

The excess energy term is described by a two-point cor- 
relations between atoms and introduces elasticity, crystalline 
symmetry and gives rise to interactions between topological 
defects within solid phases. Considering particle interactions 
truncated to second order correlations in CDFT, this term can 
be written as 
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where denotes all possible combinations of the two par- 
ticle correlations between components i and j with i,j — 

A,B,C. 

Following previous alloy PFC models (3]|6l, we define a to- 
tal mass density p = J^i Pi an d the total reference mass den- 
sity as p° = J2i Pi- Following Provatas and Majaniemi ll24ll 
and Greenwood et al. [6|, the concentration of each compo- 
nent i is defined as c,- L = p.J p and the corresponding refer- 
ence compositions are chosen as c° — p°/p°- For conve- 
nience, we define a dimensionless mass density of the form 
n = p/p° — 1. From mass conservation, we have J^. q = 1 
(or cc = 1 — ca — cb and c° c = 1 — c° A — c° B ). With these 
definitions, we re-write the free energy in terms of n and the 
{ci}. In doing so, some approximations are made for conve- 
nience. To avoid sharp density peaks, ideal free energy terms, 
expressed in n, are expanded to fourth order in the limit of 
small n. Also, we consider length scales where variations in 
concentration which are, to lowest order, slow compared to 
those of the density field, which varies on atomic scales. In 
this limit, coarse-graining of the free energy or equations of 
motion, as in Refs. [24-27], makes terms whose integrand is a 
function of slow c t fields multiplying the fast n field, vanish. 
Similarly, second order correlation terms containing combina- 
tions of the Ci can be approximated by gradient terms in Cj. 

With the above approximations at hand, the 3-component 
PFC energy functional for species A, B and C can be shown 
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where AF m ; x denotes the ideal entropy of mixing, 



(4) 



Debye-Waller-like prefactor, modulated by an effective tem- 
perature a and a transition temperature UMj lfl3l . 



Mj g 



(8) 



The k = mode of C 2 j have been omitted for simplicity but 
can be added through a constant. Its omission merely shifts 
the value of the average density and its effect can be modelled 
through the parameters of the model. 
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In Eq.|4] ?y and x are parameters introduced to control the vari- 
ation of the ideal free energy density away from the reference 
density p° . A parameter u> is introduced to correct the entropy 
of mixing away from the reference compositions c° A and c° B . 
The gradient energy coefficients a a and a b set the scale and 
energy of compositional interfaces. The parameters i], 
a a and ac have been shown to have contributions from from 
higher order correlative interactions |27|. Here, we treat these 
parameters as free coefficients to match the free energy func- 
tional quantitatively to desired materials properties. 

The correlation function, C e //, in Eq. |4] formally includes 
contributions from cross correlation functions of the form C\j 
in the excess energy. Extending the formalism of Ref. [6| 
C e f f to the case of ternary alloys, we define an effective cor- 
relation function in terms C 2 according to 
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where the coefficients Xi are polynomial functions, which in- 
terpolate between two-body correlation kernels of the pure 
species, weighing each by the local compositions. The or- 
der of the coefficient Xi vary depending on the number of 
components in the system and their order must be such as to 
smoothly interpolate from one correlation kernel to another. 
The coefficients satisfy Xi + X 2 + X 3 = 1 at all composi- 
tions. In this study they are defined by 
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X 3 = 1 - 34 + 2c\ - 3c| + 24 - 4c A c B (l -ca- c B ). 



B. Model Dynamics 

The dynamical equations of motion for each density field 
follow dissipative dynamics with stochastic noise in each field 
l28l . When re-writing these equations in terms of a total den- 
sity and two concentration fields, and neglecting the noise 
terms in each equation, we arrive at equations analogous to 
those used by Elder and co-workers in their alloy model 131 . 
namely, 
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dc B 
dt 
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M n , M CA and M CB are dimensionless mobility coefficients 
for density and compositions fields, which are, in principle, 
functions of the density, composition fields and temperature. 
In this study, at a given temperature, the mobility coefficients 
will be set to constants (i.e., equal to 1 for all mobilities). 



(7) 



The C 2 (k) are defined in Fourier space by peaks at kj, 
which corresponds to the inverse of interplanar spacings of 
the main reflection from the j th family of planes in the unit 
cell of the crystal structure favoured by component i. Each 
peak in reciprocal space is represented by the following Gaus- 
sian form of width ctj, the height of which is modulated by a 



IH. EQUILIBRIUM PROPERTIES 

We examined the equilibrium properties of a 2D Al-Cu- 
Mg system by constructing isothermal phase diagrams at two 
temperatures, i.e., a = 0.155 and a = 0.04, correspond- 
ing to solutionizing and quench/ageing temperatures, respec- 
tively. These temperatures are chosen to be well below the 
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solid-liquid coexistence (solidus) temperature. The coexis- 
tence lines at each temperature are obtained by performing 
a common tangent plane construction on the free energy sur- 
faces of solid phases expanded around the reference density 
(n = 0). The free energy surfaces of solid phases with square 
symmetry are calculated using a two-mode approximation of 
the density fields, given by, 



n t (r) 



E 2-rriki j .?/ at 

1=1 
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where the subscript i denotes a particular solid phase with a 
lattice spacing aj, and the index j counts the number of modes 
1 • • • Ni in the i-phase. Aj is the amplitude of mode j and the 
index I counts the Nj reciprocal space peaks representing the 
mode j. kij is then the reciprocal lattice vector corresponding 
to the index I in family j, normalized to a lattice spacing of 1. 

The free energy surface for each phase can be calculated 
as a function of the two composition variables, cc u and cm 9 , 
by substituting the above approximation of density field into 
Eq. Q and integrating over the unit cell. This resulting free 
energy of the crystalline phase is then minimized for the am- 
plitudes Aj. The minimization methodology is described in 
more detail in ref . [13, 23) . 
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FIG. 1. The Al-rich side of an isothermal cut (i.e., at 400°C) from 
the experimental phase diagram of Al-Cu-Mg system reprinted from 
Ref. (29 1 ; Marked by the dashed circles are the Al-rich, (Al), Cli- 
nch, (9), and Mg-rich, (/3) phase concentration regions considered 
for reconstruction by the PFC model phase diagram, as shown in 
Fig. [2] The dashed line represents the compositional boundary for 
Cu and Mg in our phase diagram calculations. 

In the Al-rich corner of the experimental Al-Cu-Mg phase 
diagram, shown in Fig. [T] there is a binary eutectic transi- 
tion between the Al-rich (Al)-fcc phase and an intermediate 
phase 9 (containing w 32.5at.% Cu) with a tetragonal crystal 
structure. The (Al)-6 system has a small solubility for Mg, 
reaching a maximum of w 2at.% Mg near the Al-rich side. 
Following the phase diagram in Fig.[T[ adding more Mg to the 
(Al)-9 system leads to the formation of a series of intermedi- 
ate phases, such as S, T and ft The latter is the cubic ftphase 
in the binary Al-Mg system with Mg-content of ps 38.bat. %. 
In our 2D model, we construct a phase diagramthat that maps 
onto the ternary eutectic system of (Al)-6-j3, with all solid 



phases having square symmetry but differing in Cu and Mg- 
content. The lattice constant (and thus the reciprocal space 
peaks) of 9 is interpolated between that of pure Al and Cu at 
32.5a£.%Cu. This calculation is performed also for the lattice 
constant of f3 considering 38.5crf.%Mg in the Al-Mg system. 
The free energy of the solid phase is generally calculated with 
a variable lattice constant weighted by concentrations cc u and 
CMg using the interpolation functions defined in Eqs. [7] The 
polynomial fitting parameters in Eq. Q (namely -q, \ and lu) 
and the width of various peaks (aj) in the correlation kernel 
are then selected so as to obtain approximately the same 
solubility limits for Cu and Mg in the solid phases as those 
in the experimental phase diagram. The parameters used are 
given in the caption of Fig. [2] 

Fig. [2] shows the free energy landscapes of solid and liq- 
uid along with the corresponding phase diagrams constructed 
for solutionizing and ageing temperature parameters of the 
model, i.e., a = 0.155 and a = 0.04, respectively. To con- 
struct the isothermal phase diagrams, the coexistence (solidus) 
lines for (Al) - 9, (Al) - ft /3 - 6 and (Al) - (3 - 6 were 
obtained by requiring that the chemical potential and grand 
potential are equal for each species in the chosen phases. For 
example, the following set of equations were solved to find 
the (Al) — 9 coexistence line: 



u {Al) = u e 



r(Al) 



(Al)JAl) _ n (At)JAl) 



c Mg C Mg 



f 6 - II 9 r 6 

J Pec, C Cu 



^c Mg C Mgi 

(13) 



where pL CQu = df/d(c Cu ) and (j, CMg = df/d(c Mg ) are the 
chemical potentials of the concentrations cq u and cm 9 , re- 
spectively. Figure |2jd) shows good agreement, in terms of 
the solubility of Mg, of the single-phase (Al) and multi-phase 
(Al) — 9 structures, in the dilute-Mg part of the experimen- 
tal phase diagram of Fig. [T] and that of the constructed one at 
a = 0.04. 



IV. CLUSTERING SIMULATIONS 

This section presents simulation results of the clustering 
phenomenon in Al-Cu-Mg system in the form of microstruc- 
tural and compositional evolution of clusters, and results are 
compared with relevant experimental evidence gathered from 
the literature. 

Using the calculated equilibrium properties discussed in the 
previous section, simulations of clustering were performed on 
a 2D rectangular mesh with grid spacing dx = 0.125 and 
time step dt — 10. The size of the grid was 4096 x 4096 
gird spacings (equivalent to 512 x 512 atoms). Each atomic 
spacing was resolved by 8 mesh spacings considering lattice 
parameter of 1 for a 2D square structure. To solve the dynam- 
ical equations, a semi-implicit algorithm was used in Fourier 
space for higher efficiency. The initial conditions were cho- 
sen to include the quenched-in dislocations in the bulk crys- 
tal, which are proposed to play a dominant role during the 
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FIG. 2. Solid and liquid energy landscapes of a square-square-square {(Al)-j3-ff) system at temperatures (a) a — 0.155 and (c) a = 0.04; 
Corresponding reconstructed phase diagrams at temperatures (b) a — 0.155 and (d) a = 0.04; The parameters for ideal free energy and 
entropy of mixing were r\ = 1.4, x = 1> w = 0.005, <f Cu — 0.333 and c° Mg — 0.333. Widths of the correlations peaks are taken an = 0.8 
and aio = v2an for all phases (the required ratio for isotropic elastic constants in a solid phase with square symmetry [13]). The peak 
positions are fc 11(Ai) = 2%, k 10 (Ai) = V^fcii(jl0> feu = (2.0822)7r, k 10e = y/2kne, hup = (1.8765)7r and ft 10 /3 = y/2kup. For 
simplicity, the effective transition temperatures Guj are set to 0.55 for all familes of planes in all phases; The concentrations co-a and cm 3 are 
rescaled considering the maximum Cu and Mg-content in the 6*-phase and /3-phase according to the experimental phase diagram, i.e., » 32.5 
and w 38.5 at.%, respectively; The concentrations on the isothermal cuts are read in Cartesian coordinates. 



early stage clustering in quenched/aged Al-Cu and Al-Cu-Mg 
alloys 031 [B] El 123 EH). Following Fallah et al. flj], ini- 
tial conditions employed a a crystal lattice of uniform com- 
position, which is distorted at the quench temperature by in- 
troducing a uniform distribution of dislocations. As with the 
previous study ifTTI . the character of dislocations is not of our 
concern in this 2D study, and they are simply defined with 
edge dislocation characteristics, with the dislocation line per- 
pendicular to the surface. While, in this study, we investigate 
the dislocation-mediated mechanisms of clustering, we ex- 
pect that the proposed mechanisms will also hold for vacancy- 
assisted clustering. 

Fig. [3ja)-(d) and (e)-(h) show the PFC simulation results 
for quench/aging of Al-l.lCu and Al-l.lCu-0.2Mg, respec- 
tively, from the solutionizing temperature of a — 0.155 to 
a = 0.04. Labelled on these images are the typical stable 
clusters "a" and "a"' in Al- 1 . 1 Cu and Al- 1 . 1 Cu-0.2Mg alloys, 



respectively, that survived the growth competition among the 
other clusters. As can be inferred from Fig. [3] the number of 
observed clusters within the unit area of simulation, at each 
time step, is much larger for Al-l.lCu-0.2Mg alloy compared 
to Al-l.lCu alloy. The zoomed-in images of the area within 
the box labelled Cluster "a"' in Fig. [^h) is shown in Fig. Qa) 
and (b) at time steps t = 1, 000 and t — 30, 000, respectively. 
Fig. |4fa) shows how initially the Cu atoms segregate into the 
areas around the dislocations. Over the simulation time, as il- 
lustrated in Fig. |4fb), the system undergoes a process of rear- 
rangement and/or annihilation of dislocations within the ma- 
trix, leading to the formation of Cluster "a"'. 

The selected alloy compositions for this study (i.e., Al- 
l.lCu and Al-l.lCu-0.2Mg) are within the single-phase re- 
gion of (Al) at the solutionizing temperature, while they 
are located within the two-phase region of (Al)-0 at the 
quench/ageing temperature (i.e., a = 0.155 and a = 0.04, 
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referring to the phase diagrams in Figj2jb) and (d), respec- 
tively). With these concentrations, during clustering simula- 
tions, we expect to observe Cu-rich clusters evolving with the 
composition of the equilibrium 0-phase (as can be seen on the 
final microstructures shown on Fig.[3jd) and (h)). Also, hav- 
ing set the square symmetry for all species, the lattice param- 
eter will be the only structural factor expected to change as 
a cluster evolves in the matrix during the 2D clustering sim- 
ulations. Intermediate non-equilibrium phases in Al-Cu-Mg 
system, which are not predicted by the experimental phase di- 
agram, were not considered in this study. 

During the simulation of ageing process in both alloys, first 
small clusters form with a slightly higher Cu and Mg-content 
than that of the matrix. As ageing progresses, some of these 
clusters shrink in size and become depleted in Cu and Mg, 
but a few become stabilized (e.g. all the clusters shown in 
Fig. [3jd) and (h), and typically those labelled "a" and "a"' 
for Al-l.lCu and Al-l.lCu-0.2Mg alloys, respectively). In 
contrast, ageing at a temperature within the single-phase (Al) 
region for both alloys, i.e., a = 0.145, results in complete 
elimination of distortion, as expected. 



raising the Mg-content in these alloys. Notably, the increase 
in number density was more pronounced for smaller cluster 
sizes, leading to a smaller average cluster size for alloys with 
higher Mg-content. Furthermore, using the PAS technique for 
these alloys in the quenched-state, a significant increase in the 
positron lifetime was recorded at higher Mg-contents, indicat- 
ing that Mg may stabilise the free volume in the matrix (i.e., 
increasing the number density of quenched-in defects), pos- 
sibly by co-clustering with Cu lfl6l [181 . Somoza et al. lfT8l 
also observed that in Al-l.74Cu-0.35Mg alloy, immediately 
after quenching, Cu-content at the positron annihilation site 
was higher than that in AI-1.74Cu alloy. They also pointed 
out that Cu-Mg complexes acted as embryos for further aggre- 
gation of Cu resulting in higher kinetics of clustering during 
ageing of the solutionaizd/quenched Al-l.74Cu-0.35Mg alloy 
compared to that in Al-1.74Cu alloy. Rapid solute aggregation 
during the early stage of ageing of solutionized/quenched Al- 
Cu-Mg alloys has been suggested through calorimetric mea- 
surements by Zahra et al. ifJTl . 



I. Effect of Mg on the evolution of microstructure 

A typical simulation result is used to study the effect of 
adding Mg to Al-l.lCu alloy on the evolution of clusters, as 
illustrated in Fig. [3] The cluster radius R, measured in terms 
of the number of atoms, is defined by radially averaging the 
radius of the concentration field bound by a threshold of con- 
centration Cp l u defined by 



2^ c Cu 



th 

c Cu 



-Cu 



-Cu 



N 



(14) 



which neglects the small concentrations of Mg. Here, N is 
the number of mesh points within the selected domain con- 
taining the cluster, and superscript 'b' denotes bulk proper- 
ties. c b Cu defines the far field concentration of Cu within the 
above selected domain. The number density of clusters is es- 
timated by normalising the number of clusters of each size 
range within the simulation domain with respect to the unit 
cell area of a 2D square crystal structure with the lattice pa- 
rameter of 2rff: c — 2.86 x 10~ 8 cm. The number density 
of clusters and their size distribution is plotted in Fig. [5] at 
t = 30,000 time steps. The data shows that Mg alloying 
in Al-Cu alloy promotes clustering through an increase in the 
number density and a reduction in the average size of the clus- 
ters, as also noted from Fig. [3] 

These simulation results are consistent with a number of ex- 
perimental observations made on the quenched/naturally-aged 
Al-Cu-Mg alloys using a combination of TEM Qg), PAS ESI 
[181 and 3D atom probe technique (APT) EUEU. With APT 
technique, Ringer et al. ||2T1 and Marceau et al. |16| found 
pre-precipitate Cu-Mg co-clusters of «3-20 atoms are dis- 
tributed within the solid solution matrix at an early-stage age- 
ing of the solutionized and quenched Al-l.lCu-(0.2-1.7)Mg 
alloys. In particular, Marceau et al. |[T6l l noted a marked in- 
crease in the number density of clusters of various sizes by 



2. Evolution of cluster composition 

The compositional evolution of cluster "a"' (marked on the 
images shown in Fig. [3jd)-(h), as it grows in the solution- 
ized/quenched Al-l.lCu-0.2Mg alloy, is illustrated in Fig. [6] 
The cluster composition is estimated by averaging within a 
circle of radius R, which is determined by a threshold Cu- 
content (i.e., Cp n as defined by Eq. 14). From the data of 



Fig. [6] we find that the Cu-content continuously rises towards 
its equilibrium value in the #-phase (i.e., c^? ~ 32.8ai.%), 
as specified by the constructed phase diagram (Pig. [2j. Mean- 
while, Mg-content and Mg/Cu ratio reach their maximum be- 
fore they continuously drop close to their equilibrium val- 
ues (i.e., c2 g ~ 0.29at.% and (^) eq ~ 0.009) in the 6- 
phase. Here, the formation process of a Cu-rich cluster with 
the equilibrium concentration can be divided into two suc- 
cessive steps; First, an increase in both Cu- and Mg-content 
within the small early clusters, and second, a reduction in 
the Mg-content while the cluster continues to attract more Cu 
atoms and grows in size until it forms a stable Cu-rich Cu-Mg 
cluster. The above evolution of composition during the two- 
stage clustering is investigated in more depth in the next sec- 
tion through a detailed analysis of the system thermodynam- 
ics, which is also supported by relevant data from the existing 
experiments. 



V. ANALYSIS OF CLUSTER FORMATION 

Formation of clusters is studied through the analyses of the 
energy in the system from the earliest small embryos to final 
stable clusters. First, the energetic mechanisms of composi- 
tional evolution is investigated for a small embryo forming 
around strained regions in the matrix (i.e., around the disloca- 
tions displacement fields) until it grows into a stable cluster. 
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t = 6,000 t = 30,000 
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FIG. 3. Time evolution of clusters in solutionized/quenched (a)-(d) Al-l.lCu and (e)-(h) Al-l.lCu-0.2Mg alloys at model temperature a = 
0.04; Typical stable clusters "a" and "a'" are labelled on the microstructures of the Al-l.lCu and Al-l.lCu-0.2Mg alloys, respectively; Each 
image represents a simulation domain with 512 x 512 atoms. 



at.% Cu: 



FIG. 4. Zoomed-in images of the area within the box labelled Cluster "a'" (in Fig.[3jh)) in the Al-l.lCu-0.2Mg alloy at (a) t = 1, 000 and (b) 
at t — 30, 000 time steps; Black dots indicate atomic positions and T-like symbols label dislocations. 



A. Compositional change from embryo to cluster 



The early small clusters, also being called embryos, show 
an increase in the Mg content and Mg/Cu ratio (as illustrated 
in Fig. |6|l up to « 14 atoms in radius. This phenomenon can 
be qualitatively explained by considering a metastable coexis- 



tence between the strain-free bulk matrix and the strained ar- 
eas near dislocations containing a small cluster (i.e., blue and 
red energy surfaces, plotted in Fig.[7|a) and (b), respectively). 
Following the methodology of Fallah et al. ifTTl . the effect of 
strain on the mean field approximation of the free energy of 
the system is evaluated in the solid-state, at a given tempera- 
ture of a = 0.04, and shown in Fig.[7jb). In this methodology, 
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FIG. 5. Number density vs. cluster size distribution in solution- 
ized/quenched Al-l.lCu and Al-l.lCu-0.2Mg alloys. 
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FIG. 6. Evolution of (a) composition and (b) Mg/Cu ratio for cluster 
"a'" shown in Fig.[3]e)-(h). R — defines the centre of the cluster. 



the magnitude of distortion caused by the displacement fields 
of dislocations near the cluster is evaluated and mapped onto 
an effective uniform strain. The strain magnitude is then used 
for calculating the free energy of the distorted structure shown 
in Fig. [7jd), which represents an area containing the early- 
stage development of a cluster (labelled "a"' in Fig.[3jf)). The 
concentration map of Cu is shown overlaid on the microstruc- 
ture in Fig. |7jd), revealing the segregation of Cu into the dis- 
torted areas near dislocations. The concentration map of Mg 
also follows the same path qualitatively. The strain value is ap- 
proximated with the following formula over the triangulated 
density peaks using the Delaunay Triangulation method: 



N tri 3 

EE 

i=l j = l 



(15) 



where N tr i is the number of triangles in the field, a is the di- 
mensionless equilibrium lattice parameter (the number of grid 
points resolving one lattice spacing, i.e., 8), is the length 
of the j th side of the i th triangle. The magnitude of strain for 
the particular structure shown on Fig. [TJd) is estimated to be 
e w 0.02. At a given temperature, the strain can be evalu- 
ated by calculating the peaks of the correlation kernel C^j in 
reciprocal space, at locations kj defined by a slight deviation 
from those of the equilibrium density peaks. The following 
equation defines the amount of strain: 



= \k-kj\/k 



(16) 



where index j denotes one family of planes in reciprocal 
space. 

For the alloy Al-l.lCu-0.2Mg, the common tangent plane 
construction illustrated in Fig. [7jc) reveals a set of composi- 
tions for a coexisting unstrained matrix and strained solid in 
the vicinity of dislocations. Particularly, it denotes higher con- 
centrations of Mg and Cu in the strained areas than those in 
the bulk unstrained matrix, while the two are in a metastable 
coexistence. Such phenomenon of segregation of solute atoms 
into the strained areas has been shown to be driven by a stress 
relaxation mechanism ifTTl [121 132| . Fallah et al. ifTTTl have 
shown that this mechanism reduces and/or even completely 
eliminates the energy barrier for formation of stable clusters 
in binary alloys. 

The above phenomenon is responsible for the first rise in 
the Mg-content illustrated in Fig. |6ja). It can be argued that 
initially both Mg and Cu atoms segregate into nearby dislo- 
cations to relax their stress fields while also forming a small 
Cu-Mg co-cluster. The cluster composition estimated from 
the above metastable coexistence (i.e., cc u ^ 6.7at.% and 
c Mg 0.7at.%) is similar to that measured for cluster "a"' 
at the peak of Mg-content (see Fig. |6|a)). The first rise in the 
Mg-content and thus in the Mg/Cu ratio is consistent with the 
atom probe data of Marceau et al. ifTBI showing an increase 
in the Mg/Cu ratio by the size of clusters up to sal 3 atoms 
in the naturally aged Al-l.lCu-(0.2-1.7)Mg alloys. Also, the 
PAS investigation of clustering in Al-l.74Cu-0.35Mg and Al- 
1.7Cu-1.3Mg by Somoza et al. (19) and Nagai et al. 11201 , 
respectively, showed that, right after quenching, a high vol- 
ume of the bulk crystal defects are associated with solute 
atoms fl9ll20l . 

In the second stage of clustering, while a randomly se- 
lected small cluster continues to grow among the other clus- 
ters, it attracts more Cu atoms causing the Mg-content to de- 
crease (see Fig. |6|a)). This occurs since the growing clus- 
ter releases more stress from the matrix by attracting more 
Cu atoms and moving towards the composition of the equi- 
librium 0-phase, which contains less Mg (i.e., ~ 0.29at.%, 
referring to the phase diagram shown in Fig j2jd)) than the 
highly-strained initial Cu-Mg co-cluster. This phenomenon 
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FIG. 7. (a) The free energy of the unstrained matrix; (b) The mean-field approximation of the free energy of a uniformly strained matrix at 
e = 0.02; (c) Metastable coexistence between the strain-free matrix (blue surface) and the strained regions (red surface) around the dislocations 
at a temperature of a = 0.04; (d) Area around Cluster labelled "a'" in Fig.Blf) at an early-stage of the transformation, t — 2, 000 time steps. 



is in accordance with the experimental evidence of the in- 
crease in the local Cu-content at the positron annihilation sites 
(i.e., vacancy-Cu-Mg complexes) during ageing of the solu- 
tionized/quenched Al-l.74Cu-0.35Mg alloy ESI [HI. More- 
over, the PAS investigation of clustering during early-stage 
ageing of Al-l.74Cu-0.35Mg fH and Al-l.7Cu-l.3Mg lESl 
alloys showed that the bulk crystal defects mingle more effec- 
tively with Cu-Mg [19, 20 1 aggregates, which are considered 
as more efficient positron traps than the single solute atoms. 

The above results and analysis suggest that although the 
crystal defects are strongly attached to the solute atoms (i.e., 
both Cu and Mg) immediately after quenching, they associate 
more effectively with Cu atoms than Mg during the initial 
growth of clusters upon ageing |fT9l . More specifically, Na- 
gai et al. EU1 report that defects are more effectively bound 
to Mg atoms rather than Cu atoms in the as-quenched state. 
However, during the early stages of the subsequent ageing, 
they observe clustering of Cu-rich Cu-Mg complexes along 
the dislocations. 



alloy, we first analyse the work of formation of a long-lived 
cluster. The work of formation for clustering is defined as 

Wh = 27ri?7 + 7ri? 2 (-A/ + AG S ) (17) 

where R is the cluster radius in terms of number of lattice 
spacings, A/ is the bulk driving force for nucleation of a 
cluster at a given concentration, AG S represents the strain en- 
ergy for a coherent nucleus and 7 is the interfacial free energy 
per unit length of the interface. We consider 2D in this work 
but we expect the same mechanisms reported below to hold 
in 3D. Particularly, our preliminary simulations in 3D have 
yilded qualitatively similar results, which will be presented in 
a sequel study. The above quantities are estimated below for 
clustering in Al-Cu-Mg alloys lying within the two-phase re- 
gion of (Al)-0. The evolving clusters are then Cu-rich, close 
in concentration and lattice parameter to those of the equilib- 
rium 9 phase, as shown on the calculated phase diagram in 
Fig-Hd). 



B. Work of formation 



1. Driving force for the formation of clusters 



In order to investigate the effect of trace addition of Mg on 
the clustering behaviour of a quenched/aged dilute Al-Cu-Mg 



The bulk driving force for the formation of a Cu-rich Cu- 
Mg co-cluster with the equilibrium concentration is defined 
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2. Strain energy created by clusters 



-A/ = f 



*4„ 



( c Mg 



(4 

- C Mg) 



4u) 



(18) 



The strain energy for a coherent nucleus is evaluated 
by 03: 



where superscripts 'b' and 'cV denote the bulk matrix and 
cluster "phase" quantities, respectively. Fig. |HJa) shows a 
low-Mg section of (Al)-0-(3 solid-phase free energy diagram. 
Using this free energy diagram, we compute the approximate 
driving force for clustering of a Cu-rich Cu-Mg co-cluster at 
two different levels of Mg-content (see Fig. |HJb)). The illus- 
tration of these calculations on Fig. [HJb) excludes the effect 

all difference in 
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on the driving force, assuming the over- 
trie Mg content between the matrix and final 
Cu-rich cluster is negligible. It is seen that Mg alloying in- 
creases the driving force for clustering of the Cu-rich clusters 
(see also Fig.[9]c)). In other words, addition of a small amount 
of Mg into the Al-l.lCu alloy decreases the solubility of Cu 
(i.e., from w 0.8 to « 0.6 at.%Cu by adding 0.2crf.%Mg, ob- 
tained from the common tangent construction in Fig. |HJb)), 
thus raising its supersaturation in the quenched state. 
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FIG. 8. (a) (Al)-9-/3 solid-phase free energy diagram in dilute Mg 
region; (b) Schematic illustration of calculation of the driving force 
for clustering of a Cu-rich Cu-Mg co-cluster at two different levels 
of Mg content. 



AG S = 2G A S 



2 K B 



K B + G A - 



(19) 



where 



S = {acu - aAi){ccu - c Cu ) + (a Mg - a A i)(c Ma - c Mg ) 

(20) 

is the misfit strain and G A and Kb are 2D shear and bulk 
moduli, respectively, calculated from PFC 2D mode approx- 
imation (6j. Following Greenwood et al. (6), in the limit of 
small deformations E6l . the free energy of different strained 
states in a 2D crystal of square symmetry is evaluated through 
substituting their respective coordinate transformations into a 
two mode approximation of the density field and integrating 
over the bounds of the strained unit cell. The elastic constants 
Cu, C12 and C44 (C12 = C44 = Cn/3) are then extracted 
through fitting the resultant free energy to parabolic expan- 
sions in displacement fields. In two dimensions, the shear 
and bulk moduli can be simply defined as G A — C44 and 
Kb = Cll + Cl2 (radially averaged for square symmetry O), 
respectively J2]- 

We calculate the bulk modulus for the Cu-rich cluster of the 
equilibrium concentration and the shear modulus for the Al- 
rich matrix for varying Mg-content, according to the approach 
presented in Ref. [33 1. Figures |9|a) and (b) show the effect of 
Mg alloying on the elastic moduli of Al-l.lCu-.TMg alloys 
and the misfit strain, 6, respectively, created by clustering of 
a Cu-rich cluster of the equilibrium concentration. For higher 
Mg contents, the combination of a reduction in the bulk and 
shear moduli of the matrix and an increase in the misfit strain 
created by the formation of a Cu-rich cluster results in higher 
strain energy values (as illustrated in Fig. [9jc)). 



3. Surface energy of clusters 

Assuming a low dislocation density in the system, the in- 
terfacial free energy is taken to be solely chemical, neglect- 
ing the structural contributions 11341 . Following Cahn and 
Hilliard ||35l , the interfacial energy between a Cu-Mg co- 
cluster and the Al-rich matrix is evaluated by the following 
analytical form which uses the composition-dependent mean- 
field free energy, f(c Cu , c Mg ),- 



acuif - f") 



OLCu dccu 



dc, 



'Cu, 



(21) 



where ac u an d aug are gradient energy coefficients for Cu 
and Mg, respectively, both set to 1 in this study. The term 
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AI-llCu cMg(at.%) Al-l.lCu cMg (at.%) 

FIG. 9. The effect of Mg on (a) the shear and bulk moduli, Ga and Kb, respectively, (b) the misfit strain 5, (c) the bulk strain energy, 
AG S , and the driving force of formation, A/, and (d) the surface energy associated with Cu-rich Cu-Mg co-clusters in an Al-rich matrix of 
Al-l.lCu-xMg alloys. 



dCMg 

dccu 



is estimated by 



8/ 

Of 
' a Mg 



and the variation of Cj^q with 



respect to ecu is approximated by a linear interpolation be- 
tween the bulk and the equilibrium cluster compositions (i.e., 



CMg - 



Fig. |9|d), increasing the amount of Mg in Al-l.lCu-xMg al- 
loys decreases the surface energy of the Cu-rich clusters. 

Fig. 10 a) depicts the evaluation of the work of formation 
in Eq. 17 which combines the effects of the computed val- 



'-(ccu — c cu)y As illustrated in 



ues of the surface and strain energy and the driving force of 
clustering. As can be inferred from this figure, the addition of 
Mg to the Al-l.lCu alloy lowers the energy barrier height and 
the critical size of nuclei for the formation of a stable Cu-rich 
Cu-Mg cluster. 



4. Effect of dislocations 

The effect of dislocations on the formation of clusters can 
be described augmenting the work of formation by the follow- 
ing form ifm 

W d = W h - AG sr + AG d 

= 2nRj + 7rR 2 (-Af + AG s )-AG sr + AG d (22) 



AG sr is the stress relaxation term due to segregation of solute 
into dislocations [36|, described here by 



AG sr = V 2 Xd,EAln(R), 



(23) 



where A = 



G A Sbf 



E 
2G A 



1, E = t? B 3 A is the 



2D Young's modulus, r\ is the linear expansion coefficient 
with respect to concentration, \d represents the change in 
the diffusion potentials due to concentration, and we intro- 
duce E6? to represent a weighted average of the magnitude 
of Burger's vectors around the dislocations accompanying the 
cluster and a is the lattice parameter. The prefactor rfxdEA 
accounts for the reduction in strain energy due to solute seg- 
regation around a dislocation [37]. Following the work of 
Cahn and Larche ll37l . which approximates the analysis of 
King et al. [38] and works well for dilute solutions, we obtain 
the following approximation for the term rj 2 Xd in an isotropic 
ternary system 



V 2 Xd 



VMgVMgg, 



d 2 f 



ZlJCuljMg-Q 



d 2 f a 2 / 



( a 2 / 

\_ dccudcM, 



1 da 



and rjMg 



1 da 



(24) 



The estimated 



where 7 1Cu - - 55 ^ — -img - "su^ 
values of the term rj 2 Xd show an increase from 0.0547 to 
0.0823 for the bulk matrix of Al-l.lCu and Al-l.lCu-0.2Mg 
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FIG. 10. Work of formation for (a) coherent/homogeneous (evaluating Eq.|22} and (b) dislocation-assisted (evaluating Eq.[l7]for one disloca- 
tion, i.e., E6f = 1) clustering of Cu-rich Cu-Mg co-clusters in an Al-rich matrix of Al-l.lCu-xMg alloys; (c) The variation of numerically 
evaluated total energy, AG to t, due to the formation of clusters "a" and "a'", within the boxes marked on the images shown in Fig. |5Ja)-(d) 
and (e)-(h), respectively. 



alloys, respectively. This implies that adding Mg to the Al- 
1. lCu alloy enhances the stress relaxation within the distorted 
matrix through segregation of solutes (both Cu and Mg) into 
areas near the dislocations. 



The term AGd in Eq.(22 1 accounts for the increase in the 
total energy of the system due to presence of dislocations. Its 
form is approximated by 



AG d = CA, 



(25) 



where £ is a prefactor of the order 10, representing the average 
amount of energy per dislocation core P9l . 

Eq.((22)> for the work of formation for dislocation-assisted 
clustering in Al-l.lCu and Al-l.lCu-0.2Mg alloys is plotted 
in Fig.[T0|b). This figure shows that for the case of a strained 
area around one dislocation (i.e., £&? = 1) the energy barrier 
height for clustering of Cu-rich Cu-Mg co-clusters of any Mg- 
content is lower than that of the homogeneous clustering (as 
depicted in Fig. [TUfa) and (b) with barrier heights labelled Bh 
and Bd for homogeneous and dislocation-mediated clustering, 
respectively. 

The total work of formation, AGtot, was also computed 
numerically by measuring the change in the grand potential 
within a box engulfing clusters "a" and "a"' during their for- 
mation and growth in the bulk matrix, i.e., 



AGtot = I <*> - I 
V Jv 



If - Hc Cu -C Cu - Hc Ug -CMg ~ fJ-n-n] 



(26) 



Here, 



sf 

Sccu 



Mcm 9 = s^- g and M» = £ are dif- 
fusion potentials of concentration and density fields, respec- 
tively, and V is the total volume. The total work of formation, 
AGtot, has contributions from the surface energy and driving 
force for formation of clusters (i.e., AGtot — AG 7 — AG V ), 



both including also the elastic effects 
Fig. [Kfc), the total 
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As can be seen in 
work of formation increases with the 



The peak in the total work of formation can be explained in 
such a way that the measuring boxes (marked on the images 
shown in Fig. [3} contain only one growing cluster, i.e., clus- 
ter "a" or "a"'. Therefore, the calculated change in the grand 
potential accounts for structural and compositional changes 
during the formation and growth of the targeted cluster. It 
should be noted that while the growth of one cluster may raise 
the local free energy, other parts of the system can undergo a 
process of annihilation and/or shrinkage of sub-critical clus- 
ters, thus resulting in a decrease of the total free energy of the 
system. 

According to our simulation data, in the presence of mul- 
tiple dislocations, typical clusters "a" or "a"' continuously 
grow until they become stable. This suggests that at each 
sub-critical cluster size (i.e., below ps 30 and s» 50 atoms 
in radius, equivalent to the radius of clusters "a" and "a"', 



respectively, at the peak of energy shown in Fig. 10 c)), the 
system is sitting at a local energy minimum. Consistent with 
the previously proposed mechanism ifTTI . locally straining a 
sub-critical cluster due to local accumulation of the magni- 
tudes of a collection of dislocations Burger's vectors creates a 
sufficiently large S6| (as illustrated in Fig.Qb)), which even- 
tually leads to complete elimination of the energy barrier. It 
is thus thermodynamically favorable for the sub-critical clus- 
ters to accumulate solute atoms from the matrix and grow in 
size. However, since the accumulation of a sufficient number 
of dislocations into the neighbourhood of a cluster is a statisti- 
cal occurrence, not all sub-critical clusters will be long-lived, 
and some of them may even reverse their growth and eventu- 
ally disappear (if not grow beyond the critical size). 

Past the critical size of the stable cluster, even the local free 
energy associated with the growing cluster starts to decrease 
due to the dominant role of the driving force over the interfa- 
cial energy. The observed peak in the numerically evaluated 
total work of formation, AG to t, is consistent with the esti- 
mated energy barriers obtained through analysis of the work 
of formation for homogeneous and dislocation-assisted clus- 
tering (see plots in Fig. 10 a) and (b)). One can easily con- 



growth of clusters to a maximum value and then decreases. 



elude that addition of Mg to the Al-1. lCu alloy reduces energy 
barrier and the critical size of the stable cluster thus leading to 
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a higher clustering rate and a finer distribution of clusters. 

VI. CONCLUSION 

In summary, we utilized a ternary extension of the alloy 
phase field crystal model of ref. |[T3l to simulate and analyze 
the atomistic mechanisms governing the early-stage cluster- 
ing phenomenon in ternary alloys. Our previous energy anal- 
ysis full of dislocation-mediated clustering in the binary Al- 
Cu system was extended to include the effects from adding a 
ternary element. Particularly, consistent with the existing ex- 
periments, we showed that Mg alloying in Al-Cu-Mg system 
refines the final microstructure. 

The detailed analysis of the system energetics for differ- 
ent levels of Mg-content in quenched/aged Al-Cu-Mg alloys 
revealed that the addition of Mg increases the effective driv- 
ing force for nucleation (i.e., A/ — AG,,), decreases the sur- 
face energy, 7, and enhances the dislocation stress relaxation, 
AG sr , associated with the Cu-rich Cu-Mg co-clusters. This 
in turn ensured a higher rate of nucleation thus leading to a 
finer distribution of clusters, a phenomenon which was also 
observed experimentally [[T6l [TBI . Furthermore, we showed 
that the simulation results for the compositional evolution of 
the early clusters are in accordance with the previously ob- 
tained chemical composition data through PAS [fT6l [T8l l20l 
and 3D APT [16| analyses. Particularly, through free energy 
analysis of the areas engulfing the typical stable clusters, we 
showed that it is thermodynamically favourable for the small 
(sub-critical) early clusters to contain more Mg than the larger 
(over-critical) Cu-rich clusters. Such phenomenon led us to 
define a two-stage clustering in Al-Cu-Mg system, compris- 



ing an initial enrichment of the cluster in both Cu and Mg fol- 
lowed by attraction of only Cu atoms leading to a reduction in 
the Mg-content, while also moving towards equilibrium con- 
centrations in both components. 

To our knowledge, this is the first ternary phase field crys- 
tal (PFC) model that elucidates the role of ternary elements 
on the phenomenon of clustering. Through a self-consistent 
coupling of the diffusive and elasto-plastic effects, our sim- 
ulations have shed light on the poorly understood atomistic 
mechanisms through which the chemical composition impacts 
on the early-stage clustering phenomena in ternary alloys, 
which are governed by crystal defects and their elastic inter- 
actions. We expect our results to hold qualitatively in 3D. 
Particularly, our preliminary 3D simulations of clustering in 
dilute Al-Cu alloys reveal dislocation-mediated formation and 
growth of early clusters, the results of which will be presented 
in a sequel study. The energy-based methodology presented in 
this study will be extended to 3D in the future to include more 
complex crystal structures and/or chemical compositions. Fi- 
nally, while we have investigated the clustering behaviour in 
the presence of dislocations, we expect that the salient fea- 
tures of the mechanisms proposed in this work will hold for 
vacancy-assisted clustering and ageing in multi -component al- 
loys. 
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